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We consider the hyperfine structure, the atomic spectrum and the decay channels of the bound 
H + fjT system (dimuonium). The annihilation lifetimes of low- lying atomic states of the system 
lie in the 10 -1 s range. The decay rates could be measured by detection of the decay products 
(high energy photons or electron-positron pairs). The hyperfine structure splitting of the dimuonic 
system and its decay rate are influenced by electronic vacuum polarization effects in the far time- 
like asymptotic region. This constitutes a previously unexplored kinematic regime. We evaluate 
next-to-leading order radiative corrections to the decay rate of low-lying atomic states. We also 
obtain order a 5 m M corrections to the hyperfine splitting of the IS and 2S levels. 



I. INTRODUCTION 



The bound system consisting of two muons (dimuonium) can be produced in heavy nuclei inelastic scattering at 
high energies and in particle decays. The decay of the neutral i]° meson into dimuonium has been investigated 
theoretically by L. Nemenov H (n° — * dimuonium + 7). The formation of dimuonium in pion-proton collisions 
(ir~ + p — ► dimuonium + n) and by photons on nuclei (7 + Z — *■ dimuonium + Z) has been discussed in [^J. For the 
direct production of dimuonium in muon-antimuon collisions, considerable experimental difficulties associated with 
slow muon beams would have to be overcome. Another possible pathway for the production of the system, which 
we do not discuss in any further detail here, is the e + e~-annihilation (near or above the /i threshold). Dimuonium, 
once produced, undergoes atomic decay (into energetically lower atomic states) and annihilation decay (into electrons 
and photons). Because the annihilation products are hard photons and relativistic electron-positron pairs, the decays 
could be investigated experimentally by established methods of particle physics. 

In this work, we devote special attention to the annihilation decay rates of low-lying atomic states. The decay rate 
of the S = 1 ortho states of dimuonium has been evaluated by J. Malenfant in pioneering investigations The 
name "dimuonium" has been proposed in Q. This work can therefore be regarded as a continuation of the earlier 
investigations on the system. We try to refine the analysis of the ortho-states, and we investigate in addition the 
S = para states. We obtain corrections to the decay rate in next-to-leading order. 

In this work, we also investigate the hyperfine structure in next-to-leading order. We further discuss briefly the 
atomic spectrum of dimuonium. The contribution of electronic vacuum polarization to the hyperfine structure and 
to the decay rate of ortho states lies in the far time-like asymptotic region. Its observation would constitute a test 
of QED in this previously unexplored kinematic regime. Due to the small length scale of the system, the hyperfine 
structure of the dimuonic atom is influenced by the effect of hadronic vacuum polarization (at the level of one part in 
10 3 ). 

Among the exotic atomic systems, some attention has recently been devoted to pionium (7r + 7r _ -system). The 
spectrum and decay channels of pionium have been discussed extensively (see e. g. |B-R)- The production of pionium 
has been observed recently M. The formation of bound atomic systems in particle decays has been observed earlier. 
Positronium has been detected in the decay 7r° — ► Ps + 7 H . The bound system of a pion and a muon (tt/j,) has been 
observed in the decay K\ — > (71-^ atom) + v JTo| |. 
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In contrast to the short-lived n + ir~ system, where c • r is in the range of 1 /Ltm, c • r is in the range of 1 mm 
for dimuonium, so that the atom could leave the beam target after production. In any of the cases discussed the 
production rate is expected to be proportional to the atomic wave function at the origin \ip n i(0)\ 2 = a 3 610/ [8 nn 3 ] 
(we use relativistic units in which % = c = 1 and a = e 2 ). Hence only low- lying states with zero orbital moment (5 
states) are expected to be produced. 

This paper is organized as follows. We first study the spectrum of the system including radiative corrections (Lamb 
shift, Section ||). We then proceed to an investigation of the hyperfine structure in Section [II. As the highest 



accuracy can presumably be obtained in decay rate measurements, we devote special attention to the decay channels 
which are investigated in detail in Sections |v| and [v|. 

II. THE ATOMIC SPECTRUM OF DIMUONIUM 

The dimuonic atom is an analogue of the hydrogen, muonic hydrogen and positronium atoms. The main feature of 
muonic hydrogen (/iH) which differs from properties of hydrogen is the order of magnitude of the Lamb shift relative 
to the fine structure splitting (see e. g. pT| ). The fs splitting in /iH and H (also in positronium and dimuonium) 
is of the order of a 4 m r , where m r denotes the reduced mass of the system, which is roughly equal to the mass of 
the lighter particle in a system with a heavy nucleus. In hydrogen the Lamb shift of the order of a 5 m r , whereas in 
muonic systems, the Lamb shift is of the order of a 3 m r , because of the large effect of vacuum polarization. The above 
considerations are parametrically true for all states, but numerical coefficients turn out to be largest for 5 states. 
Because the coefficients to the Lamb shift are negative (attractive Uehling potential) , 5 states are significantly lowered 
in energy in muonic systems. In Fig. [I] we present an overview over the spectrum of hydrogen, fiH, positronium and 

In both hydrogen and fiH, the hyperfine structure is a smaller effect of order a 4 (m 2 /M) (M denotes the mass 
of the nucleus). This is different in positronium. Positronium has a more complicated spectrum than the hydrogen 
atom. Due to the larger magnetic moment of the "central" particle (and its smaller mass), and due to virtual pair 
annihilation processes, the effects which are responsible for the fine and hyperfine structure and for the Lamb shift 
have the same order of magnitude: a 4 m r . For positronium as well as for hydrogen, the fine structure is defined as 
the separation of levels with the same principal and angular orbital momentum quantum numbers. 

In positronium, different quantum numbers have to be used for the classification of levels than in an atom with 
a heavy nucleus. For a system with a heavy nucleus, the total electron moment j — I + s is conserved within the 
central field approximation. States are classified as nljF. In positronium the total spin S = s± + S2 is an exactly 
conserved quantity. States are classifyed as nlSF (denoted as n 2S+1 lF), where 5=1 for an ortho and 5 = for para 
states. In all systems considered we denote the total angular atomic moment by F. Although the atomic spectrum 
of dimuonium is determined by a large Lamb shift of order a 3 m r (as is the case for /uH), the classification of levels 
is like in positronium (quantum numbers nlSF). 

There exists no close analogy of dimuonium with pionium, because pionium containes only spinless particles, and 
the spectrum includes neither a fine nor a hyperfine structure. Due to the strong interaction, the hadronic pionium 
system lives only a very short time, and the atomic structure cannot be investigated (because the decay rate exceeds 
the rate of atomic transitions). By contrast, for P-states in dimuonium, the atomic transition nP — > n 'S [pi < n) 



dominates over the annihilation decay rate. The lifetime of the atomic levels is considered in section ID 

The main contribution to the Lamb shift in dimuonium is caused by the one-loop electronic vacuum polarization. 
Due to the small spatial separation of the two muons in the system, the charge of the two particles, which is 

larger than the effective charge to be observed at larger distance, is less screened and causes a deviation of the energy 
levels of the order a 3 m r . For hydrogen, the vacuum polarization enters at the order of a 5 m r , but for dimuonium, 
two powers of a are compensated for by the muon-electron mass ratio (or, by employing a coordinate space picture, 
by the smaller separation of the particles). 

We now turn to a brief discussion of the atomic levels of dimuonium. The effect of electronic vacuum polarization on 
the spectrum is investigated. Any vacuum polarization insertion in the photon line can be represented as a substitution 
of the form 

- - I" ds p(s) 2 1 (1) 
q z + 1 e 7T J So q z — s + IE 

in the photon propagator. The integration has to be performed from the threshold sq of pair production of the loop 
particle. For electronic vacuum polarization, the spectral function is given by (see for instance [n2[, p. 323) 
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The substitution 



v 2 = 1 



Ami 



brings the one-loop vacuum polarization integral into the form fl3|| 

1 



o 



q + i e 



1 



1 , w 2 (l-« 2 /3) 

dv — ; o o — — — 

7r ./ 1 — v z g z — y- + ie 



(3) 



(4) 



For space-like momentum transfer, this can be Fourier-transformed into coordinate space and yields the Uehling 
potential, 



Vu(r) 

We have introduced the notation, 



1 v*(l-v*/3) 
dv : - exp (— Xr) 



1 



A 



2 m„ 



(5) 



(6) 



The correction to the energy due to the diagram in Fig. ^ is the main contribution to the Lamb shift. The evaluation 
of the matrix element of the Uehling potential on the non-relativistic wave functions leads to the results listed in 
Table [|. The Rydberg constant for the dimuonic atom is given by 



En = 



1406.6133(5) eV, 



(7) 



using the recommended values for o^ 1 = 137.0359895(61) and m M = 105.658389(91) MeV The approximate 

Lamb shift values are (neglecting higher order corrections corrections which are estimated to be suppressed by an 
additional factor of a/w): 



£(15) = -0.49 eV, £(25) = -0.058 eV, £(2P) = -0.0014 eV. 



(8) 



In contrast to muonic hydrogen (/iH), dimuonium is a purely leptonic system. Therefore its spectrum is not influenced 
by nuclear structure effects. 



III. HYPERFINE STRUCTURE OF DIMUONIUM 



The fine and hyperfine structure in dimuonium and positronium in lowest order a 4 m are given by the following 
formula (see for details textbooks ]Q^2,Qg| and articles fl7|-|l9|| ) , 



EpJnlSF) — — - — ~ m + a 4 m 
4n z 



11 , S s ,i {7 



where 



1 - S l>0 

64n 4 ' n 3 ^6 "'' u ' 4(2Z + T 



f 3 ; + 4 for F — 1 -A- 1 

(i+Tj (21+3) iU1 ^ — 1 + 



B u = < -J(ITT) for F = l > 



-i|nj for F = l-1. 
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2n 3 (2^ + 1) 



(9) 



(10) 



Here m is the mass of the constituent particle, i. e. m = m e for positronium, m — m M for dimuonium. Note that the 
above formula determines the position of energy levels in positronium up to the order of a . Radiative corrections 
in Ps enter at the order of a 5 , at the order of a 3 in dimuonium. However, for the fine and hyperfine structure 
of dimuonium in lowest order, the above formula remains valid. The results for the hyperfine structure order in 
dimuonium in leading order are 
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E^(nS) = Ep/n 3 , where E F = —a 



7 



m. 



0.175eV = 4.23- 10 8 MHz. 



(11) 



The hyperfine structure interval E^ s (nS) arises from the exchange and annihilation diagrams (see Fig. |3|) , which 
contribute 4/7 Ep and 3/7 Ep to the Fermi energy Ep, respectively. The triplet n 3 Si states are energetically higher 
than the singlet n 1 ^ states. For P states, the a 4 corrections imply the following order of the spectrum, in decreasing 
energy: n 3 P2, n x Pi, n 3 Pi, n 3 Po (we follow here the usual spectroscopic notation n 2S+1 lp). 

We now turn to the evaluation of radiative corrections to the hyperfine structure splitting for S states of order 
a/irEp (corresponding to a/nE^ s ). Some contributions are the same for positronium and dimuonium, some are 
specific to the dimuonic system. The contributions which are alike for both systems are depicted in Fig. |J These con- 
tributions originate from the anomalous magnetic moment of the electron in the transverse photon diagram (denoted 
as "(g-2)-T"), from recoil corrections (denoted by "Rec"), from the vertex correction to the annihilation diagram 
(denoted by "Vert-A"), from the muonic vacuum polarization ("VP-/1-A"), and from box-diagrams corresponding to 
two-photon annihilation (denoted by "2A"). In this paper, due to the multitude of corrections considered, we prefer 
to denote the contributions by short abbreviations, which we hope are self-explanatory, rather then single letters, 
in order to enhance the clarity of presentation. The results for the above corrections are listed in Table |] (first 5 
contributions, above the separating line). The sum of corrections to the Fermi energy of order a/nEp amounts to 

Era 



AE Ps (nS) 



a 



32 6, n 

n 2 

21 7 



7 



(12) 



where the imaginary part is entirely due to paradimuonium. It corresponds to the two-photon decay of the para 
system HIT) 



1 



r(°Vs ) 



-(0) 



2n 3 
a 5 m M 



= n 3 ■ 0.6021 • 10" 12 s. 



(13) 



Now we turn to the evaluation of corrections specific to dimuonium (the relevant diagrams are depicted in Fig. |^). 
We first evaluate the correction due to vacuum polarization insertions in the Coulomb photon line. For the dimuonic 
atom, the correction to the wave function due to electronic vacuum polarization has relative order a/ir and therefore 
modifies the hyperfine splitting in a/ir relative order. The order of magnitude of this correction is specific to the 
dimuonic atom (small length scale of the system, see Section 0. For "VPC-T" (see Fig. |J), the contribution is given 
by the matrix element 



A£vpc-t = 2x^ (V>|5(r) G(Etp) Vu\*I>), 
6 mf, 



(14) 



where 



\tpn)(lpn\ 

E^j — E n 



(15) 



is the reduced Coulomb Green function (the pole of the reference state is excluded from the sum over all intermediate 
states). Because atomic momenta are of the order of am^ in the system under consideration, it is sufficient to carry 
out the calculations with non-relativistic wave functions. The reduced Green function is a function of two coordinates 
ri and T2- Due to the appearance of the <5-like potential the Green function is needed only for r\ = 0. It can be 
expressed in closed form Q , 



Gis(£is;0,r) 



ami 2e" zi/2 



47T 



2zi(lnzi + C) + z 2 -5zi-2 



(16) 



and 



G 2 s(£ 2 s;0,r) 



,2 _-»/2 



4tt 2z 2 



4z 2 (z 2 - 2)(lnz 2 + C) + z 3 - Vizi 



6z 2 



(17) 



where C = 0.5772... is the Euler constant, z„ = (2am r r)/n and m r denotes the reduced mass of the system 
(m r = for the dimuonic atom). The radial integration can be performed by standard formulae. We finally 

obtain a one-dimensional integral over the v parameter of the Uchling potential. For the IS* ground state, we obtain 
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AE VPC -t(1S) = -E f 

7T 



, 4x 2 (l-v 2 /3) 
dv — , x 



(2 + x VT^v 2 ) 3 



7x vT 
(0.605) 



VT 



VT 



in 1 



■VT 



(18) 



where 



rn e 



For 25, a analogous evaluation yields 



A£'vpc-t(2S') 



= 1.50886. (19) 

= _ (0.523)- -f. (20) 

The analytical calculations were partially performed with the computer algebra system Mathematica p2fl . The 
contribution due to the VPC-A diagram can be evaluated by considering the relative contribution of these diagrams to 
first order hfs (4/7 and 3/7, respectively). An independent evaluation is done by utilization of a spectral decomposition 
of the reduced Coulomb Green function. If that representation of the Green function is chosen, then the calculation 
is closely related to Q. In order to simplify the calculation, we observe that the perturbation can be expressed as 
a modification of the wave function at r = 0. The discrete and continuous spectra give distinct contributions. The 
results of the respective terms (discrete and continuous spectrum) for the IS and 25* states are 



a 4 
A-Evpc-t(IS') = — Ep • 2 • - 

7T 7 



and 



AE 



VPC 



-t(2S) 



a E F 4 
— • 2 • - 

7T 8 7 



A^2s(0) 

iMo) 



= - [0.04 + 0.56] • £ F 

7T 



= - [-0.14 + 0.66] 



(21) 



(22) 



The continuous spectrum contribution (second numerical term in Eq. ( |2 1|) and (|22|) is very large. It can be understood 
in the following way. If one would try to omit the energy of the intermediate state in the denominator of the expression 



AV-(O) 



d 3 rG(£^,0,r)V' C /(r)V>(r) 



d A r 



E 



V>nt(0)^,(r) 

Eil, — En 



Vu(r) ip{r) 



in the sense of E^ — E n — > then the result should include a divergence in the sum over states. In view of the 
equation Ylnteib V'niW VC;( r ) = ^( r ) — "0(0) , 0*( T *) (completeness of the spectum), this is seen as follows: 



/ d?rY,-^^°)^iir)Vu(r)ip{r) 

nl ^ 



Vt/W- ty\Vu\il>) 
E,u 



and Vfj(O) is a diverging quantitiy. Its apperance means that virtual intermediate states with large wave numbers k > 
a are important for the convergence of the integral for A^(0), and they lead to a numerically large contribution. 
The final results for the VPC-T and VPC-A contributions are listed in Table || 

The correction due to the vacuum polarization insertion into the transverse photon line (VPT) is proportional to 
the matrix element 



M 



\V 2 Vu\i>) 



- / dv 

7T JO 



dv 



v 2 (1 - v 2 /3) 

1-v 2 
v 2 (1 - v 2 /2>) 

1-v 2 



d 3 P 



tp(r) 



V 2 

d 3 k 



exp 



(2tt) 3 J (2tt) 



■ip*(p) 



2 m e r 
s/T^v 2 
Air a (p — k) 



(p - k) 2 +4m 2 /(l - v 2 ) 



(23) 
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The matrix element M. is evaluated both in momentum space and in coordinate space with the same result. We 
obtain for the IS state 



M(1S) 



dv 



v 2 (l-v 2 /3) 2X + am p 



■In Jn " l-v 2 (X + am p ) 2 



where A is defined in Eq. (^). The final results for VPT (after numerical w-integration) are 
AE VPT {IS) = - 0.345 • E F and AE VPT {2S) = - 0.355 • E F . 



(24) 



(25) 



The results are included in Table |fi|. 

The term due to the muonic vacuum polarization in the photon line in the annihilation diagram (denoted "VP-/X- 
A" ) is known, since it has the same relative magnitude as that of electronic vacuum polarization in the positronium 
system. The electronic part (denoted as VP-e-A) is found from the well known asymptotic behavior of the vacuum 
polarization p],E2| 



AEvp cA (nS) 



a 3 

7T 7 



Ep 

n 3 



(26) 



for q 2 = (2 m^) 2 . The contribution due to the electronic vacuum polarization in the annihilation is the second largest 
correction to the dimuonium hyperfme splitting in a /it relative units (see Table |l| result for VP-e-A). The time- 
like vector q = (2m p ,0) is lying in the far time-like asymptotic region for the electronic vacuum polarization. The 
logarithmic term which appeared originally by the replacement In (q 2 /m 2 ) — > In ( — q 2 /to 2 ) (going from space like 
to time-like virtual photons) leads to the imaginary part. It can be ascribed to the decay of the ortho state into free 
electrons (/i + /i~ — > 7 — > e + e~) 



r(°)(n 3 5i 



-(()) 



6 rr 



n 3 ■ 1.806- 10 _12 s. 



(27) 



In contrast to positronium, the parastate and orthostate lifetime have the same order of magnitude in dimuonium 
(a 5 rrifj). 

The hadronic contribution (denoted by VP-h-A) to the vacuum polarization is found as the sum of four terms. The 
main contribution results from a pionic loop. We follow here the approach in p5| . The spectral function is of the 
form 



P(s) 



4 m 2 ) 



2 \3/2 



12 s 5 / 2 

where the pionic form factor is used in the form given by Gounaris and Sakurai [p6[ 

N 



D 1 +D 2 -iD 3 



The quantities A, D\, D2 and D3 are given by 



N = m p + drrip T p , 



(28) 



(29) 



(30) 



3 m 2 ^ m p + 2 k p m p m 2 m p 



7r k 2 2 m T 2 7T k p tt k% 



0.48, 



(31) 



D x = m 2 - s, L> 3 = nip T p (k(s)/k p \ m p /^/s, 



k(s) 2 (h( S )-h p ) + k 2 p h\m 2 p ) (m 2 p - S ) 



(32) 
(33) 



where h' denotes the derivative of h, and the functions k and h arc defined as 
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H«) = ;V S -'1< h(s) = -=■ In 

2 7r ys \ im^ 



(34) 



with the special values k p = k{m 2 p ), h p = h{m 2 p ),T p = 150.7(1.2) MeV, m p = 768.5(6) MeV |§. We give results for 
the 15 state only in the sequel. The 1/n 3 scaling is easily restored in the final result. The contribution from the 
pionic vacuum polarization is given by 



A£, + ,-(15) = -~ 

7T i 



Ami / ds P{S) 



' Ural ^ m l- s 



Ep = — (-0.055) • E F . (35) 

7T 



In the simple p-meson pole approximation, where p(s) — At: 2 / f p S(s — m 2 p ) with f 2 /(An) — 2.2 p7[ , we have 

A^ OT+W - (15) w AS P (15) = - (-0.050) £ F . (36) 



This result agrees with the full pion form factor of ( |29| ) to about 10 %, so it is justified to consider mesonic resonances 
of higher energy in the pole approximation only. The higher energy mesonic contributions are due to the w and 
(j) resonances. These resonances are not included in the Gounaris-Sakurai form factor and are treated separately. 
Estimating the coupling constants as f 2 /(Ai:) = 18(2), /J/(47r) = 11(2) []25|,p7[ and given the meson masses of 



>4> 

rn u = 782 MeV, m = 1019 MeV |15j, the results are 

OL OL 

AE U (1S) = -(-0.006) -Ep and AEJlS) = - (-0.005) • E F . (37) 

The background above 1 GeV is estimated by assuming a form of 

, . R , _ <j(e + er — > hadrons) . . 

p(s) = —, where R= \ ^—^ , (38) 

3 s cr(e + e — > /i + fi ) 

for the spectral function with a branching ratio R m 2 (constant) below y^i = 4 GeV and R w 4 above -y/s = 4 GeV 
(see p. 190). Integrating from an estimated lower threshold of s t h ~ (lGeV) 2 we obtain 

OL 

A£> (15) w - (-0.014) • £ F . (39) 

7T 

Summing all contributions and restoring the 1/n 3 scaling, we have as the contribution from the hadronic vacuum 
polarization 

OL Ep 

A£Vp-h-A(n5) = AE V + AE^ + AE^ + AE> = - [-0.080(9)] • -4 ■ (40) 

We estimate model-dependent uncertainties in the hadronic vacuum polarization to be of the order of 11%. The 
hadronic term is included in Table || (result for VP-h-A) . 

The sum of all corrections (see Table ||) to the hyperfine splitting in dimuonium amounts to 

OL OL Ep 

A£ h fs(15) = -0.689(9) • E F and A£ hfs (25) = - 0.556(9) • — ^ . (41) 

7T 7T 8 

In the final results for hyperfine splitting, we estimate higher order corrections to enter at the 5% level of the next- 
to-leading order contributions. We obtain 

E Ms (lS) = 4.23283(35) • 10 7 MHz (42) 

and 

E Ms {2S) = 5.28941(34) • 10 6 MHz . (43) 

Results for the Lamb shift and the hyperfine structure of low-lying levels are presented in Tables | and ||. The 
largest deviation from the scaling appeared in VPC-T and VPC-A contributions (16%), because the wave function 
(see Eqs. ( ^l| ) and (|2^)) is more sensitive to the behaviour of the potential about the origin. 

An important point for a possible investigation of the spectrum is the lifetime of levels. As one can see from Eqs. 
( |l3| ) and ( p7| ) the annihilation lifetimes of 5-levels are much shorter that the lifetime of the free muon, t p = 2.20 10~ 6 s. 
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By contrast, the annihilation decay rate for P-states includes an extra a 2 , and the lifetime is of the same order of 
magnitude as that of the free muon. 

The decay rates of excited states have to be compared with the atomic transition rates, which are also of the order 
of a 5 (see e. g. @). We obtain 

t(2P ^ IS) = 1.54- l(r n s. (44) 

Annihilation of P states is suppressed by two orders of a compared to S states due to the behaviour of the P wave 
function near the origin. Therefore, the annihilation lifetime of the 2P state can be estimated to be of the order of 
10 _7 s, and it is seen that the atomic transition, not the annihilation, determines the lifetime of the excited 2P state 
in dimuonium. 

This situation is different for S states, where the annihilation dominates over atomic transitions. The annihilation 
lifetimes of of 15 and 2S states lie between 0.6 • 10~ 12 s and 14 • 10 -12 s, and we expect that the decay rates could be 
measured via detection of the decay products. Radiative corrections to the decay rates are considered in the following 
Section. 



IV. DECAY CHANNELS OF ORTHODIMUONIUM 



The leading-order contributions to the orthodimuonium and paradimuonium decay rate can be extracted as the 
imaginary parts of the energy corrections to the hyperfine structure, 



a 5 m^j 
6 n 3 



and r (0) (n 1 S* ) = 



2n 3 



(45) 



The calculations at leading order were presented above (Eqs. (|13|, |27j) ) . The above results can also be found in [|3 
We begin the consideration of radiative corrections with orthodimuonium, and provide results for the IS* state here. 

Results for the 25 state (and IS) are summarized in Table III. The diagrams contributing in next-to-leading order 

for ortho states are depicted in Fig. [j| 

The VPC-A correction can be interpreted as a modification of the wave function at the origin. The energy shift and 

decay rate are both proportional to |i/>(0)| 2 . The annihilation diagram contributes 3/7 Ep to the first-order result for 

the hyperfine splitting. Hence, the VPC-A diagram yields a correction of 



Ar VPC -A(l 3 5 1 ) 



a 

TT 
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■ 0.454 



r (o)(i35 1 )--i.06-r(°)(i 3 5 1 ) 



(46) 



to the decay rate. Analogous considerations are true for the Vert-A correction, because the diagrams consist of 
separated blocks, and hence the correction to hfs and to the decay rate can be traced back to the same matrix element 
of the wave function. We obtain 



ArVert-A^Sl) 



.r(°>(n 3 £i). 



(47) 



The diagrams VP-fi-A, VP-e-A and VP-h-A have to be interpreted as modifications of the photon propagator. 
Because the energy shift is proportional to the amplitude of the propagator, but the decay rate is proportional to its 
square, we have to multiply the relative correction to hfs by a factor of 2 in order to obtain the relative correction to 
the decay rate. Hence, the VP-/Z-A, VP-e-A and VP-h-A diagrams yield a total correction of 



a 

TT 



16 








~~9 ' 


6 \ m e ) 



10 



- 0.37(4) 



• r(°VSi) 



to the decay rate. 

We now consider the one-loop radiative corrections to the electron line (ReA) and emission of a photon by an 
electron, i. e. bremsstrahlung (BeA). The sum of BeA+ReA can be easily obtained from the diagrams in Fig. ^ 
for the two-loop electron polarization correction to the hyperfine splitting, which is completely determined by the 
asymptotic behaviour of the two- loop vacuum polarization (see e. g. |l3(]). For the correction to the hyperfine splitting, 
we obtain 



Am 2 and ((3) = 1.202. 



is the Riemann £ function of argument 3. This correction is of relative order 



where q 2 

a 2 1 tv 2 with respect to Ep and was therefore not considered in Section III. The imaginary part of this contribution is 



just the result for the sum of the diagrams ReA+BeA. The correction to the decay rate can now be evaluated easily. 
The sum of ReA+BeA yields a correction of 



3 a 

4 7T 



r(°VSi) 



relative to the leading order result (cf. Eqs. (26 j45|)). 

The last term in relative order a/n is the three-photon annihilation. The result is known from orthopositronium 
calculations Q (see also 



2a b 1 7T 2 



A^A^Sr) = o 



9 



m u = — 



3 v 



9) 



• r(°)(n 3 ^). 



(49) 



The final result for the IS decay (see also Table HI) is 



Ar(l 3 5i) = - 



221 



4^ 
3 m e 36 



— + 0.68(4) 



4 (tt 2 - 9) 



r<°>(i 3 Si) 



I |(l.90) + (0.68(4)) 
- 3.74(4) -rW(l 3 5i). 



1.16)|r (0) (l 3 5i) 



(50) 



The first term corresponds to the sum of the analytically evaluated contributions Vert-A, VP-/J-A, VP-e-A and 
BcA+ReA, and has the numerical value 1.90. The second term originates from the numerically evaluated contri- 
butions VPC-A and VP-h-A. The last term is associated to the three photon decay. The result in Eq. (50) has 
to be compared with the earlier analysis of the decay rate of heavy leptoniu m B . The final result for the ana- 
lytically evaluated contributions, which is identical to the first term in Eq. (pOF, whose analytical expression is 
(4/3) ln(2m M /m e ) — 221/36, is found in agreement with the pioneering investigations by J. Malenfant For the 
VPC-A correction, Malenfants results are in slight numerical disagreement with ours (see Eqs. (77,78) in [Q). We 
presume this disagreement can be traced to the fact that Malenfant has calculated the VPC-A correction with free 
Green functions, whereas the evaluation in this work is done using bound Green functions. In this context it is 
important to note that in the limit of am^/2m e — > 0, our result is agreement with that of Malenfant. This can be 
seen as follows. The VPC-A correction may be rewritten as a correction Atp(0) to the wave function at the origin. 
After this reinterpretation, we find 



A ^(0) = - ( ^n + 0( K 2 )) m for k=^ 
7r \ 16 / Zm e 



0. 



(51) 



which is in agreement with Eq. (80) of Q. We can therefore conclude that in the limit of weak binding (k — » 0), 
Malenfants result is in agreement with ours. However, we hold the view that bound Green functions should rather be 
used for the VPC-A correction. The atomic momentum in dimuonium is of the order of am^/2 = nm e (n ps 0.75). 
This momentum is close to the mass of the loop particles (electrons and positrons) of electronic vacuum polarization. 
These particles determine the radius of the Uehling potential. Therefore, some of the momentum integration for 
the VPC-A correction is performed in an area about m e , where the bound Coulomb Green function cannot be 
approximated by the free Green function (because the effect of the binding Coulomb potential, in momentum space, 
is inversely proportional to the square of the momentum transfer). This consideration should explain the slight 
numerical disagreement for the VPC-A correction between this work and the result in B. 



The sum of next -to-leading order corrections to the decay rate for the 2S state is (see Table III) 



Ar(2 3 Si) = - 3.60(4) • r^(2 3 5i) . 



(52) 



We estimate the higher order corrections to be suppressed by an additional factor of a compared to the next-to- 
leading order result. We obtain the following results for the decay rate of orthodimuonium, taking into account also 
the uncertainty from our model of the hadronic vacuum polarization: 



r(l 3 S'i) = 1.79073(77) ■ 10~ 12 s and r(2 3 S' 1 ) = 14.3305(59) ■ 10~ i2 s 



(53) 







V. DECAY CHANNELS OF PARADIMUONIUM 



The diagrams contribu ting to paradimuonium decay in next-to-leading order a 6 m are presented in Fig |9| results 
are summarized in Table [V. The vertex correction term (Vert-2A) is equivalent to the corresponding correction for 
parapositronium pOf, 



.2 



Ar V c r t-2A(n 1 5 ) = -- E-JEll^VSo) • (54) 

TT 4 

The correction to the wave function caused by the VPC-A diagram modifies the decay rate of the para system in the 
same way as the ortho system (see Table IV), and introduces a deviation from the 1/n 3 scaling. 

For the para state, there exists another correction to the decay rate corresponding to the decay into a photon 
and an electron-positron pair. The result for the A2e correction can be obtained in the following way. We consider 
the one-loop vacuum polarization insertion into a 2-photon annihilation diagram (Fig. [l0|). Because the vacuum 
polarization insertion, evaluated for a real photon, must be equal to zero (gauge invariance of vacuum polarization) , 
the imaginary part of the diagrams in Fig. [l0| leads to the result we need. We can use Eq. (||) for the parametric 
form of the vacuum polarization insertion, and we consider the s-integration as the final one. The integrand is now 
equivalent to the imaginary part of the 2A diagram for the hfs, but with one of the photons having a finite mass s 
(s corresponds to the sum of the four-momenta of the emerging electron-positron pair) . The vacuum polarization 
insertion normally fixes a gauge (the Landau gauge) for the virtual photon, because the polarization insertion is 
proportional to the transverse projector. However, as it was demonstrated in plL it is possible to substitute in Eq. 
(^J) any covariant gauge. For convenience, we choose the Feynman gauge. As a result we have the expression 

Ar A2o (n 1 5 ) - - [ dsp(s)T^( S ,0), (55) 



where r'°'(s, 0) is the decay rate to one real photon and a virtual photon with mass s. In order to obtain the 
correction in relative units, we divide by the first-order result, which is given by r(°)(0, 0). The correction relative to 
the first-order result r'^fW-So) is given by 



r^ifn^o) tt J ' K ' r(°)(o,o) ' 

An additional factor 2 appears because the insertion of a vacuum polarization operator doubles the number of non- 
equivalent diagrams contributing to the imaginary part. For the logarithmic coefficient we may neglect s in r(°)(s,0) 
and approximate r(°)(s, 0) — > r(°)(0, 0), and use the asymptotic form of the spectral function (cf. Eq (^)), 

p(s) — > — for s — ► oo . 
3 s 

We can thus easily obtain the logarithmic coefficient, 

Ar A2c( ^ )_ f A2H2m ^ ds ± = 2 i ln ai^i ln ^ 



r^VSo) Urn.)* 3 s 3 ml 3 



The full result requires a more detailed analysis of the dependence of r(°)(s, 0) on s. It differs from the approximate 
analysis presented above only by an additive constant. The final result of the calculation is 

AT A2e ^S ) = ^(Un 2 -^- 1 4) A-r ( °V*o). (57) 



tt \3 m e 9 

An independent evaluation of the A2e-correction using the standard S-matrix formalism is used to veri fy t he result in 
Eq. ( |57| ) . Treatment of the Dirac currents involved in simplified by application of a symbolic program |3^] developed 
for high energy physics calculations by A. Hsieh and E. Yehudai. Some care must be taken during evaluation, because 
one cannot assume the electrons as massless in the final states (the result else includes a logarithmic divergence in 
the electron mass). Proper regularization of the relevant expression then leads to the result in Eq. 



The final result for the next -to-leading order radiative corrections to paradimuonium decay (see Table IV) is 



Ar(l 1 5 ) = -4.79- r<°> (l^o) . (58) 
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For the 2S state, we have 



Ar(2 1 5 ) = - 4.65 .r<°> (2 1 S Q ) 



(59) 



Estimating higher order corrections to enter at the level of 5% of the next-to-leading order contributions, we obtain 
the following theoretical values for the paradimuonium decay: 



t(1 1 S'o) = 0.59547(33) • 1(T 12 s and r(2 1 S , ) = 4.7653(25) • 1CT 12 s . 



(60) 



VI. CONCLUSIONS 



We evaluate next-to-leading order corrections to the spectrum, to the hyperfme splitting and to the decay rate of 
low-lying levels of the dimuonic system. The results for the spectrum are given in Section |l[ We observe that 
for 2P states, the atomic decay into the IS state dominates over annihilation processes. This would facilitate 
experimental observation of the atomic transition, if dimuonium atoms can be produced in quantities sufficient to 
carry out spectroscopic measurements. We evaluate the hyperfme splitting of IS and 2S states in next-to-leading 
order. The results are 



and 



Ems(1S) = 



Ehf s (2S) = 



1 



1 



-0.689(9) 

7T 



-0.556(9) 

7T 



~2 a m » 



12 



4 „, 

■ a m. 



4.23283(35) • 10 7 MHz 



5.28941(34) • 10 6 MHz. 



(61) 



(62) 



We present a complete evaluation of all next -to-leading order radiative corrections to the_lifetime of both the ortho- 
and para-state of the dimuonic atom. In leading order, we reproduce the known results 



6 71 d 



and r(n 1 S' ) 



2n a 



as imaginary contributions to the hyperfme splitting (Eqs. ( |l3| , p7|)). 
The results in next-to-leading order are for orthodimuonium, 



r(l 3 S' 1 ) = [l + - 3.74(4)1 r(°)(l s 5i) and T(2 3 S' 1 ) 

7T 



1 



a 

7T 



3.60(4)1 r^(2 3 S*!), 



(63) 



(64) 



where the primary theoretical uncertainty is due to hadronic vacuum polarization. The lifetimes of orthostates are 
given by 



r(l 3 5i) = 1.79073(77) ■ 10~ 12 s and t{2 3 S 1 ) = 14.3305(59) • 10" 12 s . 
For paradimuonium, we obtain 

r(i 1 5 ) = 



1 + -4.79 

7T 



T^^So) and T(2 1 5 ) = 



a 

1 + - 4.65 

7T 



r(°)(2 1 s , o). 



The lifetimes are given by 



r(l 1 S'o) = 0.59547(33)- 10~ 12 s and r(2 1 5 ) = 4.7653(25) • 10~ 12 s . 



(65) 



(66) 



(67) 



We estimate higher order QED corrections to enter at the 5% level of the next-to leading order corrections considered 
in this work. 

Lifetimes in the 10~ 12 s range can be measured by established methods of particle physics via detection of the decay 
products (electron-positron pairs in case of orthodimuonium and two photons in the case of paradimuonium). We 
stress that accurate decay rate measurements can be accomplished with fewer individual atoms than would be needed 
for spectroscopic investigations. 

One of the ways to investigate the hyperfme structure of 15* or 25* states could be based on the observation 
of the interference between paradimuonium and orthodimuonium in an infrared frequency field at resonance (the 
radiofrequency field would mix the two states and thus yield a modified decay rate of the statistical sample). We 
conclude that the dimuonic system offers the possibility to observe quantum electrodynamic effects in a previously 
unexplored kinematic region. 
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state C C 

IS -0.15 -0.49 eV 

2S -0.018 -0.058 eV 

2P -.000043 -0.0014 eV 



TABLE I. Contributions to the Lamb shift due to the electronic vacuum polarization. C is given in relative units, 
AE — 2 C Eo, where Eq is the Rydberg constant for the dimuonic atom. 
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diagram 



C(1S) 



8C(2S) 



(g-2)-T 

Rec 
Vert-A 
VP-^-A 

2A 



0.571 
-0.857 
-1.714 
-0.381 
0.263 



0.571 
-0.857 
-1.714 
-0.381 
0.263 



VPC-T 

VPT 
VPC-A 
VP-e-A 
VP-h-A 



0.605 
0.345 
0.454 
1.483 
-0.080(9) 



0.523 
0.355 
0.393 
1.483 
-0.080(9) 



Sum 



0.689(9) 



0.556(9) 



TABLE II. Corrections of relative order a/% to the dimuonium hyperfine splitting. All contributions are explained in 
the text. The corrections (g-2)-T, Rec, Vert-A, VP-fi-A and 2A contribute to the positronium hyperfine splitting just as 
in dimuonium. The remaining contributions (VPC-T, VPT, VPC-A, VP-e-A, VP-h-A) are specific to the dimuonic system. 
Results are given in relative units AE = a/% C Ef- 
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diagram 



C(1S) 



8C(25) 



Vert-A -4.00 -4.00 

VPC-A 1.06 0.92 

VP-^-A -1.78 -1.78 

VP-e-A 6.92 6.92 

VP-h-A -0.37(4) -0.37(4) 

BcA+ReA 0.75 0.75 

3A 1.16 1.16 



Sum 3.74(4) 3.60(4) 



TABLE III. Corrections of relative order a/ir to the orthodimuonium decay. Results are given in relative units 
AT{n 3 Si) = a/7rCT (0) (l 3 Si). 
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diagram 


C(1S) 


8C(25) 


Vert-2A 


-2.53 


-2.53 


VPC-2A 


1.06 


0.92 


A2e 


6.26 


6.26 


Sum 


4.79 


4.65 



TABLE IV. Corrections of relative order a/ir to the paradimuonium decay. Results are given in relative units 
ArVSo) =a/irCT ( - \l 1 So). 
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Hydrogen 



2 2 P 



3/2 



2 2 S 1 / 2 



2 2 Pi/ 2 



2 2 S l 



2 2 P 



3/2 



2 2 R 



1/2 



/2 



Positronium 



2 1 P 1 

2 3 P 1 



2 3 P 



2 1 P 1 
2 s P 1 



2 3 S 1 
2 1 S 



FIG. 1. Overview over the spectrum of n — 2 levels in hydrogen, muonic hydrogen (pH), positronium and dimuonium 
(fi + The double lines denote the hyperfine structure splitting of the levels in hydrogen and muonic hydrogen. Dashed lines 
denote S = 0-states, full lines denote S = 1-states. 



FIG. 2. (Electronic) vacuum polarization insertion in the Coulomb photon (main contribution to the Lamb shift in muonic 
systems). The dashed line denotes a Coulomb photon. Bold fermionic lines denote muons, thin lines denote electrons and 
positrons. 
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FIG. 3. Transverse photon exchange diagram and (time-like) photon annihilation diagram. Both of the diagrams contribute 
to the hyperfine structure of S states in leading order. The zig-zag line denotes a transverse photon, the wavy line denotes the 
full photon propagator. 
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^> <^+^> <c - L0 - 

v n , J 

Vert-A 




VP-n-A 




FIG. 4. Known corrections up to the order a/% Ef to the hyperfine splitting of dileptonic systems (positronium, dimuonium) . 
Diagrams are explained in the text. For the dimuonic atom, the corrections depicted here contribute to hfs, but there are 
additional terms specific to the dimuonic system which need to be taken into account. The subtraction of lower order (L.O.) 
contributions is necessary for some of the diagrams in order to prevent double counting. The bold fermionic lines denote muons. 
The direction of time in all diagrams is from left to right. 
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VPC-A 





VP-e-A 



VP-h-A 



FIG. 5. Corrections of order a/ir Ef to the hyperfine splitting specific to the dimuonic atom. Diagrams are further explained 
in the text. 
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FIG. 6. Decay channels of orthodimuonium (one-photon decay into an electron positron pair, OM) and paradimuonium 
(two photon decay, PM). The decay channels depicted yield the main contribution to the decay of the system (r of order 
a 5 m M ). Bold fermionic lines denote muons, thin lines denote electrons and positrons. 
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L.O. 



Vert-A 




VPC-A 



' ^ 

VP-u-A 

>6< 

VP-e-A 



VP-h-A 





VPC-A 



ReA 



BeA 



+ 5 sym. diag. 



3A 



FIG. 7. Order a/ir corrections to the decay channels of orthodimuonium. The 5 symmetrical diagrams originate from the 
symmetrization of photon wave functions. 
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FIG. 8. Evaluation of bremsstrahlung and electron vertex corrections to orthodimuonium decay as imaginary part of the 
two-loop vacuum polarization insertion in the photon line. 
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Vert-2A 




FIG. 9. Order a/ir corrections to the decay channels of paradimuonium. The symmetrical diagrams originate from the 
symmetrization of the photon wave functions. 
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FIG. 10. Evaluation ol corrections to the paradimuonium decay caused by the production ol an electron positron pair (A2e 
correction) . 
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